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I. INTRODUCTION 
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Theoretical studies on heavy baryons containing a single heavy quark have gained pace recently as a result of the 
experimental progresses in the last few years (see [l[ and references therein). The presence of a heavy b or c quark 
i in these baryons makes them more attractive, theoretically. They can provide information on the structure of the 
QCD and its parameters as well as clues on new physics effects. For this reason, it is important to understand their 
properties as precisely as possible. Couplings of these baryons to pions are one of the important properties, also 
because their pion cloud can significantly modify their properties. 
. . , In light of the experimental developments, the mass spectrum of these baryons has been studied extensively via 
QCD sum rules method both in the finite mass limit (see e.g. Q) and heavy quark limit Q, in the heavy quark 
effective theory Q and different phenomenological models (see for example (5l-ll0l|V Electromagnetic [Lll4T7j and strong 
[TlUm interactions are two of other important characteristic of hadrons that can be used to probe their structures. 

In the present work, the coupling of the Sq baryons, Q = c or 6, to pion is calculated in the light cone QCD sum 
rules (LCQSR) framework. This problem has been studied using chiral perturbation theory in |2lj . Note that, the 
same framework has also been used to calculate the strong coupling constants of the light mesons with heavy baryons 
p^| - [20j . The layout of the paper is as follows: in the next section, we obtain sum rules for the corresponding coupling 
constants in the framework of LCQSR. Section 3 contains our numerical analysis and discussion. 



II. LIGHT CONE QCD SUM RULES FOR THE ps Q s Q ^ COUPLING CONSTANT 

To study the coupling constant in LCQSR method, the following correlation function is studied: 

n = i J | T{rtetf{x)fmp(p)} | o), (i) 

where r\^nii (x) is the interpolating current of El 1 * baryon containing the quarks q\qiQ {q% = u, or d) and T stands 
for the time ordered product. The correlator can be calculated in terms of hadronic parameters inserting a complete 
set of hadronic states. This expression is called the phenomenological or physical side. It can also be calculated in 
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theoretical side in terms of QCD parameters via operator product expansion (OPE) in deep Euclidean region, where 
— p 2 — > oo and — (p + q) 2 — > oo. QCD sum rules for the considered coupling constant is obtained matching both sides 
of the correlation function and applying Borel transformation in order to suppress contribution of the higher states 
and continuum. 

First, we focus our attention to calculate the physical side of the correlation function. The phenomenological side 
of the sum rules for the heavy baryons is similar to the light baryons [22T ] . For this aim, we insert a complete set of 
hadronic state having the same quantum numbers of the interpolating currents into the correlator. As a result, we 
obtain 



jj (0 I m Q I s Q ( P2 )) (go(gi) I vz Q I 0) 

n = ^ s (T lQ (p 2 )TT(q) I £q(pi)) 2 T 1 



P% - ml 



Pi ~ m T.r 
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where pi = p + q and p 2 = P an d ■■• denotes the contribution of the higher states and continuum. The matrix element 
creating baryon from the vacuum is parameterized as: 



(0 I V^ Q I £q(P: s)> = As Q us Q 0, s), 



(3) 



where As Q is the residue and us Q is the spinor for the Sq baryon. The remaining matrix element appearing in Eq. ^) 
defines the coupling constant g£ Q £ Q7r through 



(E Q (p 2 )7r(g) I £q(pi)) = 9^Q^ Q nu{p 2 )ij5u{pi). 



Using the Eqs. Q and ((4]) in Eq. ([2]) we obtain the phenomenological side of the correlator as: 
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In principle, one can choose any of the structures, ^ ^75 and ^75, existing in Eq. ([5]). Our calculations show that the 
structure p 1 $75 leads to a more reliable result. 

After obtaining the physical side, we proceed to calculate the QCD side of the correlation function. For this purpose, 
we use the interpolating currents in the following general form: 
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riyud — —=e 



a be 



u F a CQ b Wd c + p\u T CfQ b \d c 



lCd b bV-/3 QlCj 5 d b u. 



a be 



ui CQ b 7 V + [u A a CYQ b )u c -[Q 1 a Cu b 7 V - P[Q a C^u b u. 



(6) 



In Eq. ([5]) the Q denotes the heavy quarks b or c, /3 is an arbitrary parameter with /3 = — 1 corresponding to the Ioffc 
current, C is the charge conjugation operator and a, b and c represent the color indices. 

To proceed in our calculations in QCD side, we need also the propagators of the heavy and light quarks. They are 
given as (23j : 
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16tt 2 x 2 ^ y ' ^ -~^y—'< 4lr 2 x 2_ 
The free light and heavy quark propagators in Eq. ([7} are given in x representation as 
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where i*Q are the Bessel functions. Besides the propagators, the matrix elements of the form {it{q)\q{xi)Tiq(x2)\fy 
are also required. Here r, represents any member of the Dirac basis i.e., {1, j a , a a p/^/2, *7s7a, 7s}- The matrix 
elements (?r (q) I q (x-\ )Tj q (x-i ) \ 0) are parameterized in terms of the pion light cone distribution amplitudes and they 
are given explicitly as [2J, [25[ 



(7r(p)|^)7 M75 g(0)|0) = -iUppj due^ x [^{u) + —mlx 2 K{u) 



PX Jo 



(Tr(p)\q(x)i l5 q(0)\0) = ^ / duj Upx ip P {u) 



(Tr(p)\q(x)a a pj 5 q(0)\0) = -/j, n (l - p,^) ( Pa xp - P px a ) j due lupx ip a {u), 



(Tr(p)\q(x)a^~/ 5 g s G a p(vx)q(0)\0) = 
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PlSPn [ gua ~ ~ < <P" Xa +Pa X v) 



+ PpPu Ulna ~ ^{PfiXa +p a X ll ) S j 
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(TT{p)\q(x)^ig s G a0 (vx)q(O)\O) = Pfl (p a Xp - P pX a )±-f w ml J Vae^+^^V^a,) 



+ 



P/3 [ g^a - —(PftXa + p a Xy) 



where ^ = f n - 



x / Vae l(a ^ +va <' )px V ± {a l ), (9) 
_ m u +m d ^ jy a _ d a -d aq( i as §{\ — o.q — otq — a g ) and the ip^[u), A(it), B(u), tpp(u), ip a {u), 



m-a+m d ' 

T(oti), A±(ai), A\\(ai), V±(oii) and Vii (cKi) are functions of definite twist and their expressions will be given in the 
numerical analysis section. Using these inputs, the QCD side of the correlator can be calculated in a straightforward 
fashion. 

Having calculated the correlation function both in physical and QCD sides, we equate the coefficients of the selected 
structure from both representations and apply Borel transformation with respect to variables p 2 and (p+q) 2 to suppress 
the contribution of higher states and continuum. As a result of these procedures, the QCD sum rules for coupling 
constant gs: Q s: Qn is obtained as: 



Si) 



Q 



r 'm SQ |A SQ | .gs Q E Q 7r = / e ^ p(s)ds + e 7m*T, 



(10) 
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where the functions, p(s) and T are given as: 
1 
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(13) 



and /i(ct!i) = V||(ai), / 2 (a l ) = Vi_(aj), /s^i) = (ck*), A^i) = TCa^), / 5 (o!j) = vT(ai) are the pion distribution 
amplitudes. Note that, in the above equations, the Borel parameter M 2 is defined as M 2 — .^V^a and uq = 



Mr 



Since the masses of the initial and final baryons are the same, we can set M 2 = M| and Uo = 4. 

As it is clear from the expression of the coupling constant, we need also the residue, As Q , which is calculated using 
mass sum rules 1201 : 



, 2 -mi. IM 



e"M^pi(s)ds + e « 2 Ti 



(14) 



with 
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—^[5 + 0(2 + 5/3)][12^io - 6^ 2 o + 2</>30 - 4^ 4 i + V>42 - 12Zn(-V)], 



(15) 
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< dd >< uu > 



2M 4 AM 2 



(16) 



Note that, since only the square of the residue appears in Eq. (fl~4|) . only the magnitude (not the sign) of As Q can 
be determined from mass sum rules. This indeterminacy is also carried into the coupling constant calculation, and 
hence sum rules determine the coupling constant upto a sign. 



III. NUMERICAL ANALYSIS 



This section is devoted to the numerical analysis of the coupling constant gs Q s Q 7r- The numerical values of 
the required input parameters are given as: (qq)(l GeV) = — (0.243) 3 GeV 3 , nib = 4.7 GeV, m c = 1.27 GeV, 
m Sb = 5.805 GeV ,m Ec = 2.455 GeV, mg(l GeV) = (0.8 ±0.2) GeV 2 [27|, U = 0.131 [MUzI and TO ff = 0.135 GeV. 
In order to calculate the coupling constant, the 7r-meson wave functions are needed and explicit forms of the these 
wave functions are represented as [24j, |25[ 
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where C*(x) are the Gegenbauer polynomials 

aio 

"10 

hoi 

h w 
9o 

!)2 
.94 

The constants in the Eqs. (JXTj) and (fT8|) were calculated at the renormalization scale fi = 1 GeV 2 using QCD sum 
rules [24-26] and are given as a\ = 0, a£ = 0.44, % = 0.015, r? 4 = 10, w 3 = -3 and w 4 = 0.2. 

The sum rules for the coupling constants contain three auxiliary parameters, the Borel mass parameter M 2 , the 
continuum threshold sq and the general parameter j3. These are not physical parameters, hence the sum rules for 
the coupling constants should be independent of them. Therefore, we should look for the working regions for these 
parameters such that the dependency on these parameters are weak. At too large values of the Borel parameter 
M 2 , the suppression of the contribution of higher states and continuum is reduced, increasing the error due to the 
quark-hadron duality approximation. The lower limit of the M 2 can be obtained by requiring that the highest 
twist contribution, which is suppressed by a higher power of M 2 , should be small. Therefore, the higher states and 
continuum contributions should comprise a small percentage of the total dispersion integral. The continuum threshold 
so is not completely arbitrary and is related to the energy of the first exited state. From the numerical analysis of 
the sum rules for the coupling constants, the continuum thresholds sq are obtained as 38 GeV 2 < s„ < 42 GeV 2 and 
9 GeV 2 < Sf) < 11 GeV 2 for the E^E;,^ and E c E c 7r vertexes, respectively. In order to acquire the working region 
for ft parameter, we plot the <7£ Q £ Q7r as a function of cos# in the interval —1 < cos# < 1, which corresponds to 
— oo < f3 < oo, where tan# = (3. From the Figs. (PJ and which show the dependency of coupling constants on 
the cos#, it is seen that for —0.5 < cos (9 < 0.4 the dependence of the coupling constants on cos 6 is weak. We present 
also the dependency of the coupling constants, gs c s c Ti and .gs b s bir , on the Borel parameter M 2 in Figs. ((3]) and (j4j, 
respectively. These Figures reveal that the coupling constants are quite stable in the chosen Borel mass region. From 
these Figures, the numerical values of the coupling constants, (?£ c £ c7r and gT. b T. b -m can be extracted as follows: 

5ScSo7r = -8.0 ± 1.7, 9^^ = -H-0 ±2.1. (19) 

The errors appearing in our predictions are due to the uncertainties in the input as well as the auxiliary parameters. 

At the end of this section, let us compare our results with the existing prediction on the £(>„.+ coupling constant 
[22j for the light E baryon. In light case the s quark and in our case (baryons with a single heavy quark) the c 
or b quark is the spectator quark. These quarks do not participate in the considered stron g in teractions, so one 
expects these coupling constants (heavy and light cases) have values close to each other. In J22|, it was estimated 
that, 5S+S07T+ = —9 ± 2. The same transition for the light E baryon was also studied in [28l | and [29j]. In [28| . 
the ratio, (fe+£ 7r+ / 'snn-tt — 0.8 was obtained, which leads to the value of the coupling constant, g-£.+Y. Q -K+ — H-9 
with the experimental value of gNN-n = 14.9 [30j | . The value of the same coupling constant was given in [2{| as 
5s+s°7r+ — —11-9 ± 0.4. Comparison of our results with that of [22|, [28|, [2{| supports the expectation that the 
spectator quark does not effect the pionic coupling, significantly. 



1 

3 

21 9 

~8~ 

21 

~8~ 

7 :! 

4 
7 

4 
1, 

18 20 
1 + — al + 60773 + —rji, 

-^-6«. (18) 



-TJ4W4, 



I 774 ^ 4 20 ° 2 ' 



7 



IV. ACKNOWLEDGMENT 



The work has been supported in part by the European Union (HadronPhysics2 project " Study of strongly interacting 
matter"). 



[1] R. Mizuk, arXiv: 0712.0310 (hep-ex); M. Kreps, ECONF C 070512 (2007) 19. 

[2] F. S. Navarra, M. Nielsen, Phys. Lett. B 443 (1998) 285; D. W. Wang, M. Q. Huang, Phys. Rev. D 68 (2003) 034019; Z. 

G. Wang, Eur. Phys. J. C 54 (2008) 231; F. O. Duraes, M. Nielsen, Phys. Lett. B 658 (2007) 40; X. Liu, H. X. Chen, Y. 

R. Liu, A. Hosaka, S. L. Zhu, Phys. Rev. D 77 (2008) 014031. 
[3] E. V. Shuryak, Nucl. Phys. B 198 (1982) 83. 

[4] A. G. Grozin, O. I. Yakovlev, Phys. Lett. B 285 (1992) 254; Y. B. Dai, C. S. Huang, C. Liu, C. D. Lu, Phys. Lett. B 371 
(1996) 99; D. W. Wang, M. Q. Huang, C. Z. Li, Phys. Rev. D 65 (2002) 094036. 

[5] N. Mathur, R. Lewis , R. M. Woloshyn, Phys. Rev. D 66 (2002) 014502. 

[6] D. Ebert, R. N. Faustov, V. O. Galkin, Phys. Rev. D 72 (2005) 034026. 

[7] M. Karliner, H. J. Lipkin, arXiv: 0307343 (hep-ph), 0611306 (hep-ph). 

[8] M. Karliner, B. Kereu-Zura, H. J. Lipkin, J. L. Rosner, arXiv: 0706.2163 (hep-ph). 

[9] J. L. Rosner, Phys. Rev. D 75 (2007) 013009. 
[10] M. Karliner, H. J. Lipkin, Phys. Lett. B 575 (2003) 249. 
[11] B. Julia-Diaz, D. O. Riska, Nucl. Phys. A 739 (2004) 69. 
[12] S. Scholl, H. Weigel, Nucl. Phys. A 735 (2004) 163. 
[13] A. Faessler et. al, Phys. Rev. D 73 (2006) 094013. 
[14] B. Patel, A. K. Rai, P. C. Vinodkumar, arXiv: 0803.0221 (hep-ph). 
[15] C. S. An, Nucl. Phys. A 797 (2007) 131, Erratum-ibid; A 801 (2008) 82. 
[16] T. M. Aliev, A. Ozpineci, M. Savci, Phys. Rev. D 65 (2002) 096004. 
[17] T. M. Aliev, A. Ozpineci, M. Savci, Phys. Rev. D 65 (2002) 056008. 
[18] Peng-Zhi Huang, Hua-Xing Chen, Shi-Lin Zhu, Phys. Rev. D 80 (2009) 094007. 
[19] Zhi-Gang Wang, Eur. Phys. J. A44, 2010, 105; Phys. Rev. D 81 (2010) 036002. 
[20] K. Azizi, M. Bayar, A. Ozpineci, Phys. Rev. D 79 (2009) 056002. 
[21] Hai-Yang Cheng, Chun-Khiang Chua , Phys. Rev. D 75 (2007) 014006. 
[22] T. M. Aliev, A. Ozpineci, S. B. Yakovlev, V. Zamiralov, Phys. Rev. D 74 (2006) 116001. 
[23] I. I. Balitsky, V. M. Braun, Nucl. Phys. B 311 (1989) 541. 

[24] P. Ball, JHEP 01 (1999) 010; PBall, V. Braun, A. Lenz, JHEP 0605 (2006) 004. 
[25] P. Ball, R. Zwicky, Phys. Rev. D 71 (2005) 014015. 

[26] V. M. Braun and I. E. Filyanov, Z. Physik C 44 (1989) 157, V. M. Braun and I. E. Filyanov, Z. Physik C 48 (1990) 239, 
A. R. Zhitnitsky, I. R. Zhitnitsky and V. L. Chernyak, Sov. J. Nucl. Phys. 41 (1985) 284 [Yad. Fiz. 41 (1985) 445], V. A. 
Novikov et al., Nucl. Phys. B 237 (1984) 525. 

[27] V. M. Belyaev, B. L. Ioffe, JETP 56 (1982) 493. 

[28] G. Erkol, M. Oka, T. T. Takahashi, Phys. Lett. B 686 (2010) 36. 

[29] H. Kim, T. Doi, M. Oka and S.H. Lee, Nucl. Phys. A 678 (2000) 295. 

[30] R. A. Arndt, Z.-J. Li, L. D. Roper, R. L. Workman, Phys. Rev. Lett. 65 (1990) 157. 



8 




-14- 



FIG. 1: The dependence of the coupling constant gs c s c w on cos# for the value of Borel parameter M 2 = 7 GeV 2 and the 
different values of continuum threshold so- 
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FIG. 2: The dependence of the coupling constant (/S(,E b ir on cos# for the value of Borel parameter M 2 — 30 GeV 2 and the 
different values of continuum threshold so- 
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FIG. 3: The dependence of the coupling constant gf£ e £ C 7r on the Borel parameter M 2 for the value of arbitrary parameter 
/3 = — 3 and the different values of the continuum threshold so- 
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FIG. 4: The dependence of the coupling constant gs^E^ir on the Borel parameter M 2 for the value of arbitrary parameter 
/3 = — 3 and the different values of the continuum threshold sq. 



